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In the last few years ARCADE 2, combined with older experiments, has detected an additional 
radio background, measured as a temperature and ranging in frequency from 22 MHz to 10 GHz, 
not accounted for by known radio sources and the cosmic microwave background. One type of 
source which has not been considered in the radio background is that of fast transients (those with 
event times much less than the observing time). We present a simple estimate, and a more detailed 
calculation, for the contribution of radio transients to the diffuse background. As a timely example, 
we estimate the contribution from the recently-discovered fast radio bursts (FRBs). Although their 
contribution is likely 6 or 7 orders of magnitude too small (though there are large uncertainties in 
FRB parameters) to account for the ARCADE 2 excess, our development is general and so can be 
applied to any fast transient sources, discovered or yet to be discovered. We estimate parameter 
values necessary for transient sources to noticeably contribute to the radio background. 


I. INTRODUCTION 

Radio astronomy and astrophysics have an influential 
history, intertwined with astronomy, cosmology, particle, 
and nuclear physics. Today there are mysterious new 
sources being discovered in the radio universe. In this 
work we will focus on two recent and exciting develop¬ 
ments, namely an observed and unexplained excess in the 
radio background (ERB) and the discovery of fast radio 
bursts (FRBs). We use the latter fast radio transient 
(FRT) sources as the example for our modeling. 

The ARCADE 2 (Absolute Radiometer for Cosmol¬ 
ogy, Astrophysics and Diffuse Emission) experiment re¬ 
ported [1] an excess in temperature, after subtracting the 
cosmic microwave background (CMB) and accounting for 
other known backgrounds, in measurements at frequen¬ 
cies from ~ 3-10 GHz. When combined with previous 
data and analysis in the literature there is an excess that 
extends from 22 MHz to nearly 10 GHz. This is well fit 
by a power law in absolute temperature versus frequency 
with the following form [1] 

/ j/ S -2.5990 ± 0.0366 

r= (24.1 ±2.1 K) (- . (1) 

^ '^VSlOMHz/ 

While ARCADE 2 measures a temperature for the ra¬ 
dio background, transient sources like FRBs are mea¬ 
sured by other experiments in units of spectral irradi- 
ance or spectral flux density: Janskys,^ where IJy = 
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^ For our fellow particle physicists, one might be tempted to iden¬ 
tify the energy implicit in Watts with that in the Hertz in the 
denominator, but these energies do not “cancel.” One should 
think of the Watts as the power that is received by the radio 
telescope at a particular frequency in Hertz (and per unit area). 
The bandwidth (frequency range of observation), measured in 
Hertz as well as the observing area, will directly affect the power 



FIG. 1. A log-log plot of eq. (2) versus x, with T — 10. 


10“2® Wm“^ Hz“^. To convert these temperature mea¬ 
surements into a flux density, we show the standard dis¬ 
tribution function (i.e. Planck radiation law) for a mass¬ 
less Bose particle (here a photon) as a function of the 
dimensionless variable x = hv/ksT in Fig. 1. The en¬ 
ergy density per frequency distribution, dp/dz/, is given 
by 

dp _ 2 _ 2/U/3 1 

dv TT — 1 _ I 

with the first form written in natural units {h = c = 
kg = 1) with overall units of energy^, and the second 
is the typical textbook form with the constants restored 
and units proportional to Janskys. Note that when T is 
allowed to vary with frequency, as in the ERB data, then 
the spectral form will differ from the fixed distribution 
shown in Figure 1. 


received, and thus it is more sensible to express the radio unit as 
power per unit frequency and area (an energy flux density): the 
Jansky. 
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This excess background has persisted and has been 
difficult to understand despite different analyses (see, 
e.g., [2, 3]), and remains of unknown origin (for a recent 
study see, e.g., [4] and references therein). If the back¬ 
ground comes from additional extragalactic sources, they 
must be very numerous and faint, e.g. [2, 5]. Additionally, 
the background appears to be extremely smooth, which 
puts further constraints on the size and cosmological dis¬ 
tribution of such background sources [6]. There have 
been several different proposals for an alternative origin 
for this radio background, including dark matter [7], but 
the question is still open. 

More recently there has been the discovery, first with 
the Parkes telescope [8-10] and then at Arecibo [11], of 
a new class of radio objects, dubbed “fast radio bursts” 
(FRBs) due to their only known activity occurring on 
the millisecond Earth-time scale (FRB discoveries have 
continued as well [12]). Based upon their dispersion 
measures (DMs),^ these FRBs are believed to be extra- 
galactic in origin (although this has been challenged, see, 
e.g., [13]). The physical process powering an FRB is 
so far unknown, with proposals ranging from neturon 
stars [14], axions [15], and black holes [9, 16] to sev¬ 
eral others [17]. It is further inferred that these un¬ 
known objects release ^ 10^®-10‘*° erg during their mil¬ 
lisecond duration, corresponding to a typical peak flux 
density of > 0.4 Jy. While there are only a few ob¬ 
served FRBs (11 at this moment; an up to date list 
with measured parameters and references is available at 
http;//astro.phys.wvu.edu/FRBs/), a rough estimate 
of the rate for these events gives of order 10,000 per day 
in the whole sky out to z = 0.75 (recent studies, such 
as [18], put the rate a bit lower, though still consistent 
with earlier results and of negligible effect for our present 
analysis). 

Motivated by FRBs, we will consider general fast tran¬ 
sient radio sources as contributing to the ERB. By “fast” 
we mean that the timescale of such an event is much less 
than the time of observations (which is of order an hour 
for ARCADE 2). Such sources have not been considered 
in the past for calculating the diffuse radio background; 
we will adjust the steady state radio source calculation 
to accommodate the transient contribution. 

This paper is organized as follows. First, in Sec¬ 
tion II A, we will introduce the basic calculation for radio 
source counts for the radio background, and modify it to 
account for fast transients. While the FRB contribu¬ 
tion to the excess radio background is several orders of 
magnitude too small to make detectable contribution, a 


^ Without candidate source galaxies for FRBs thus far, the only 
way to estimate their distance is through their dispersion mea¬ 
sure. The DM is determined by the electron density along the 
line of sight, which delays lower frequency signals compared to 
higher frequencies. Given an estimate for the electron contribu¬ 
tion from the galaxy, the large DMs of FRBs lead to the inference 
of extragalactic sources. 


more detailed calculation for transient sources, e.g. sim¬ 
ilar to the calculation for supernova neutrinos, may be 
useful. In Section IIB we present a more complete anal¬ 
ysis. We then compare a general fast source to the excess 
background to derive approximate parameters for such a 
source to contribute to the excess at the percent level or 
more in Section III. Finally, in Sections IV and V we dis¬ 
cuss our results and future directions before concluding. 


II. A RADIO BACKGROUND FROM FAST 
TRANSIENTS 

Fast transient sources have not been considered as a 
contribution to the observed ERB. One might ask if 
they can contribute appreciably to the radio background. 
In this Section we will first use a simplified analysis 
to estimate the contribution from fast transients before 
providing a more comprehensive calculation. Fast tran¬ 
sient sources need to be treated slightly differently in the 
usual radio background calculation in order to properly 
account for their event rate and individual event time. 
For a recent review of FRT searches and related bounds, 
see [19] and references therein. 


A. An Estimate 

For a typical source (static, or slowly changing, com¬ 
pared to typical observation times), the contribution to 
the cosmic radio background (CRB) (see e.g. [2]) is 

r°° dN 

Bcrb= / S—dS, (3) 

‘S'min 

where S is the flux density at Earth 
(energy/(time-area-frequency), e.g. in units of Jan- 
skys), and ^ is the differential number of sources at 
that flux density. The unit for radio observations, the 
Jansky, implicitly has a time such that multiplying by 
an observation time, area, and frequency bandwidth, 
yields the total energy observed. 

For transient sources, however, there are now two 
timescales which are introduced: the duration of the 
source emission and the inverse rate of sources appear¬ 
ing in the sky. For fast sources (where “fast” means short 
compared to the observation time, e.g. ARCADE 2 here), 
the entire energy of each event is observed. Instead of a 
flux density per source, we should use an energy density: 
we multiply S by St, the time of the source mission. We 
still want the same units at the end, so we need an in¬ 
verse time, which comes from the number per flux density 
being replaced by a rate per flux density (i.e. a number 
per energy density), dN/dS dt, with units of per second 
after integration over S. In other words, the source flux 
density is to be multiplied by the duration of the tran¬ 
sient source, while the number density becomes a rate. 
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We can write this as 

The result can properly be interpreted as the incoming 
energy flux density, just as in Eq. (3).^ 

Let us consider observations at one particular fre¬ 
quency, ~ 1.4 GHz, for which there is observation of 
both the ERB^ and FRBs. The radio background can 
be expressed as lO^Jy/sr (e.g., see [2]). The estimated 
event rate of FRBs of 10^ per day in the whole sky (for 
S' > 1 Jy) corresponds to 10“^ s“^ sr“^ in units more use¬ 
ful for us. As a lower estimate on the FRB contribution 
to the ERB at 1.4 GHz, we simply multiply this rate by 
the typical flux density times the millisecond duration, 
IJy X 10“^s. Thus the FRB contribution has a lower 
bound of ^ 10“® times the observed radio background; 
this is a lower bound because the true contribution will 
be larger, from inclusion of further/dimmer events, as 
well as modeling the distribution of FRBs in redshift. 
(We will see in Section H B 1 that our more detailed cal¬ 
culation, using a model for the FRB spectrum and dis¬ 
tribution, increases the FRB background by a factor of 
~ 100.) 

To verify this very rough estimate, we will assume a 
form for the differential number per flux density, and in¬ 
tegrate backwards in time (outwards in space). (Since it 
is estimated from the DMs that the observed FRBs have 
red-shifts < 1, the red-shifting of energies, and dilation 
of time, are small.) With the assumptions that FRBs are 
standard candles observed at extragalactic distances (as 
is so far consistent with all available data), the nearby 
(z < 1 or S' > 1 Jy) distribution approximately follows 
that of a static Euclidean universe, ^ We 

normalize this distribution to the estimated ERB rate 
above, ^ dS = 10“^s“^sr“^. Inputting this as¬ 
sumption allows integration of eq. (3), and we find the 
FRB contribution matches our rough estimate of 10“® of 
the observed radio excess. 


B. A More Detailed Calculation 

To calculate the radio background due to FRTs in the 
universe, we follow the same basic procedure as in, e.g., 
the neutrino background from supernovas (see [20] for 
a review). The two main ingredients are an event rate 
(number per volume per time) evolving in redshift, i?(z). 


® Implicitly this is over the entire sky, but one can divide by 47r or 
else use a rate which is per steradian. 

^ For now we will ignore any errors in the measurements, though 
the ERB data at 1.4 GHz is consistent (within errors) with zero 
after additional source subtractions. One could use a different 
frequency for the ERB, with an assumption of the FRB spec¬ 
trum, and would find similar results. 


and the spectrum (energy per frequency) of an individ¬ 
ual event, F(y)^ with v the frequency observed at Earth. 
The incoming isotropic flux background, in units of en- 
ergy/(time • area • steradian • frequency), or Jy/sr, is 
then calculated by the following integral. 


B{y) 


^ i?(z)F(zz') 

Att Hq J g -|- z)^ Ha 


( 5 ) 


where c/Hq is the Hubble length, z^ax is the redshift 
back to which we are integrating sources, F = iy(l + z) is 
the blueshifted frequency, and flm and Ha are the energy 
densities of matter and dark energy, respectively. 

To perform this calculation we will need to model our 
radio sources. First, let us consider how the sources are 
distributed in redshift. We will employ two models, one 
where the source distribution simply scales with the co¬ 
moving volume as a standard candle (FRBs were mod¬ 
eled like this in [21]), and another which tracks the star 
formation rate. These models are shown in Figure 2. 

In the first model we assume the sources have a con¬ 
stant number density in comoving volume, then the rate 
simply scales with the fraction of comoving volume out 
to tf{z). With 




dz'- 


1 


Hq 


i(t) Hq Jq z')^ Ha 

the comoving distance, the source rate is 

D{zf 


pFHT/ _ tuFRT 

-n-CovG! — -'''Cov 


D{z,f' 


( 6 ) 

( 7 ) 


where ^ normalization to be fixed 

by the source model or observation. 

Alternatively, we model the source rate to be propor¬ 
tional to the star formation rate. There are several differ¬ 
ent parameterizations and we follow the one formulated 
in Table H of [22], but with some parameters slightly 
rounded as in [20], 


r(l + z)3-28, z<l 

7?star(^) OC < ai(U- z)“°-^®, 1 < z < 4.5 , (8) 

[a2(l-bz)-7-8, 4.5 <z 

with oi = 11.6 and 02 = 4.45 x 10® fixed by continuity. 
The transient radio source rate is then 

i?™T(z)=iVs™i?star(^), (9) 

with a normalization denoted by 

Now we turn to modeling the spectral properties of 
the fast sources themselves. As the simplest model we 
will take the spectrum of the sources to follow a power 
law (in frequency) over some range of frequency. The 
free parameters will be the spectral index, overall nor¬ 
malization, and frequency range. Since we are studying 
the radio background over observations much longer than 









4 


the source timescale, we will not need to directly assume 
this timescale. Instead, it is implicit in our normaliza¬ 
tion of the total energy. For simplicity we will assume 
the sources are either standard candles or consider all 
quantities to be typical or average values. 

We write the spectrum as 

Fil^) = Fnorm , (10) 

with spectral index a, reference frequency and nor¬ 
malization Enorm with units of energy per unit frequency. 
The normalization is fixed by integrating the spectrum 
over its frequency range, t'min to t'max, to determine the 
total energy released in an event, 

/‘^'max 

Etot.i= / di^F{iy). ( 11 ) 

^'min 

1. Fast Radio Burst Modeling 

Before considering a general (as yet unknown) source 
to fit to the excess radio background, we return to FRBs 
as an example of applying the above basic modeling. 
With the current data on FRBs, the consensus on the 
approximate rate of FRBs is about 10,000 per day in 
the whole sky, out to a redshift 0.75. We will use this 
observed rate, expressed in more useful units as 

Rfkb{z < 0.75) = 2.40 x 10"'^^ s'^ m-^, (12) 

as the integrated rate to redshift 0.75 to fix the normal¬ 
ization. 

If we assume FRBs are standard candles with a con¬ 
stant number density in comoving volume (as in [21]), 
their rate is modeled by eq. (7) with the normaliza¬ 
tion, = 4.68i?FRB(-2 < 0.75) determined by set¬ 

ting z* = 0.75 in eq. (7) and integrating out to z*. In 
our alternative model based on the star formation rate, 
eq. (8), the normalization, obtained by direct integration, 
is = 0.43Rfrb(2 < 0.75). 

The second set of parameters for FRBs is their typical 
energy output and spectrum. From the observed events 
thus far, the energy released is about 10^^ to 10^^ J, or 
10^® to 10^° erg. We use the FRB model in [21] as a 
reference point: a spectral index of —1.4 from 10 MHz to 
10 GHz, with a bolometric luminosity of 8 x 10"^'^ erg/s, 
or a total energy of 8 x 10®'* J (using an event time of 
1 ms). 

With these parameters we can calculate the FRB con¬ 
tribution to the radio background, again at 1.4 GHz. 
We find that the FRB contribution is approximately 
3 X lO”*” of the observed excess with a constant comoving 
volume event distribution, and 0.5 x 10“*" when using an 
event rate proportional to the star formation rate. 

Given the large uncertainties in what is known about 
FRBs, we can vary the parameters to some degree. From 
eq. (5) we know that any changes to the overall rate or 


FRB Rate (N 



FIG. 2. A model for the FRB rate scaled with comoving 
volume (orange, dashed) and tracking the star formation rate 
(blue, solid). The rates are normalized such that their integral 
to z = 0.75 matches the approximate observed rate. 

energy of the radio transients will proportionally change 
the calculated background contribution. For instance, if 
the rate is underestimated by a factor of a few and the 
total energy is underestimated by an order of magnitude 
(both reasonable possibilities, given the observed values 
and error bars), the FRB background contribution could 
be 10“® of the excess. While further tuning of parameters 
or modeling might improve this slightly, it seems unlikely 
that, without a drastic modification of what is known 
about FRBs, could they make up a significant proportion 
of the ERB. 


III. FITTING A FAST TRANSIENT SOURCE 
TO THE EXCESS BACKGROUND 

We will now fit the flux density of a new source of 
FRTs to the observed ERB from ARCADE 2 and pre¬ 
vious experiments. These experiments report a temper¬ 
ature beyond the contribution from the CMB tempera¬ 
ture of 2.725 K and the galactic foreground. We use “Fit 
1” of [23] to estimate the contribution from unresolved 
extragalactic discrete sources (see also [24] for resolving 
such sources). Using the Planck radiation law, Eq. (2), 
we convert these temperatures to an expected spectral 
flux density, in units of Jy/sr.® The resulting data is 
shown in Table I. We note that the data point at 1.4 GHz 
is consistent with zero given the errors on the measured 
temperature after the above subtractions. We do not in¬ 
clude data at frequencies larger than ^ 8 GHz as they are 
also consistent with zero. 

Fitting to these 8 data points, from 22 MHz to 8 GHz, 
we find a good fit to the observed excess with the contri¬ 
bution from FRTs integrated back to redshift 6. 

First, we consider the case when we model the tran¬ 
sient rate as the comoving volume. With the parameters 


® Note that the reported errors in temperature become asymmetric 
in radiance; we use the smaller of the max/min errors for each 
data point as a common error for fitting purposes. 









5 


Frequency 

Temperature (K) 

Final Temperature (K) 

Flux Density (kjy/sr) 

22 MHz 

20 355 ±5181 

13 479 

200.430 ± 77.043 

45 MHz 

3864 ± 502 

2866 

178.300 ±31.200 

408 MHz 

13.42 ±3.52 

8.107 

41.410 ± 18.000 

1.42 GHz 

3.271 ±0.526 

0.4568 

26.24 ± 32.60 

3.20 GHz 

2.787 ±0.010 

0.052 05 

2.667 ± 1.380 

3.41 GHz 

2.770 ±0.008 

0.036 62 

0.6777 ±0.4850 

7.98 GHz 

2.761 ±0.013 

0.03516 

(13.92 ± 13.90) X 10“® 

8.33 GHz 

2.743 ±0.015 

0.01725 

(73.23 ± 73.20) x 10“® 


TABLE I. The ERB data used in this analysis. The data is from ARCADE 2 and older experiments (see Table 4 of Ref. [1] 
and references therein), with the first temperature being the reported temperature before subtractions, but accounting for the 
galactic foreground. The final temperature is after the CMB, and extragalactic sources have been subtracted (the error bars 
are the same as in the previous column). The flux density is derived from this final temperature and is the value used for 
fitting. 


Jy/sr 



FIG. 3. The cosmic radio background due to fast radio tran¬ 
sients (blue, solid) with cosmological rate based on the co¬ 
moving volume, and the observed radio excess data points 
(orange) with error bars derived from the temperature mea¬ 
surements. The dashed orange line is from the continuous best 
fit function (power law in temperature) as reported by AR¬ 
CADE 2. The excess may continue up to about 10 GHz, but 
significance at these high frequencies drops and the best fit 
from the fast transients has no power at these high frequen¬ 
cies. Note that at 1.4 GHz the error bar is consistent with 
zero, where we’ve subtracted known radio sources as in [24]. 


for the simple model (the total luminosity £x rate iipRT, 
the spectrum endpoints, and the spectral index), we find 
a per degree of freedom of about 4.03/4. The spec¬ 
tral index is 0.56 and the spectrum goes from 10 MHz to 
7.3 GHz with a normalization of 4.6 x 10^^ JHz“^. Note 
that the fitted spectrum cuts off at high frequencies be¬ 
fore the final data points at 7.98 and 8.33 GHz reported 
by ARCADE 2 to have an excess temperature, but the 
reported excess of these points are highly suppressed (of 
lower significance, and consistent with zero) compared 
to the other points. Plots of the fitted model and radio 
excess are shown in Figure 3. 

For the cosmological fast transient rate modeled af¬ 
ter the star formation rate, the results are similar: The 
fitted spectrum again has spectral index 0.56, normal¬ 
ization 4.2 X 10^® JHz“^, and an overall x^ P^r degree 


Jy/sr 



FIG. 4. The same as in Figure 3 but with the fast radio 
transient rate proportional to the star formation rate. 

of freedom of 3.96/4. The htted lower frequency of the 
spectrum is the same as the previous case, but the up¬ 
per endpoint is 4.5 GHz. This model is plotted in Fig¬ 
ure 4. We note that the largest difference between this 
background model and the previous one is in the overall 
normalization of the spectrum, which needs to be about 
an order of magnitude larger in this case. 


IV. DISCUSSION 

We turn now to interpreting the parameters of the best 
fit found above. Compared to the model spectral in¬ 
dex of —1.4 for FRBs, the spectral index of the best fit 
is very flat in order to match the observed radio back¬ 
ground. And, as we know from the FRB results, we need 
an overall factor of about 7 orders of magnitude larger 
to completely saturate the radio background excess. Al¬ 
ternatively, this factor needs to be at least 5 orders of 
magnitude larger than FRBs in order to be at the per¬ 
cent level of the radio excess. 

This normalization is the combination of two main fac¬ 
tors: the event rate, and overall energy. One could also 
subdivide this into further details, such as possible re- 
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currence rate, luminosity dependence on redshift (i.e. a 
distribution of observed source energies), and so on. 

While there have been searches for FRTs in the past 
(see the recent review in [19]), much of parameter space is 
unexplored and unconstrained. Searches are also usually 
done in narrow frequency bands, and are not sensitive to 
general combinations of parameters (e.g. minimum peak 
flux density versus event time). This makes it particu¬ 
larly difficult to draw general conclusions and constraints 
on what FRT parameters are currently allowed, though 
the flat fitted spectrum here may already be constrained. 
Due to the recent discovery and study of FRBs, this set 
of parameters (bright millisecond sources at 1.4 GHz) is 
better studied. From our analysis above, it looks rather 
unlikely that FRBs, or sources very similar to FRBs, play 
a role in the ERB. 

On the other hand, there is still much more to be ex¬ 
plored. Monochromatic sources will be very easy to miss 
without searching at exactly the right frequency, for in¬ 
stance. While a single such source could not explain the 
entirety of the ERB, perhaps multiple types of sources 
play some part. It is not difficult to find parameters for 
FRTs that avoid current constraints (which, we stress, 
are typically not very general), especially rare very bright 
sources (perhaps peaking at unexplored frequencies), or 
very dim sources. A comprehensive analysis of the avail¬ 
able parameter space, and potential FRT models, is out¬ 
side the scope of this current work, but is a promising 
future direction. With instruments such as the Square 
Kilometer Array aiming to search for radio transients, 
the coming years may yield many new discoveries. 


V. CONCLUSIONS 

In this work our primary aim has been to extend the 
analysis of sources which can contribute to a radio back¬ 
ground to include FRTs. These have not been considered 
before in relation to the cosmic radio background. With 
large swaths of parameter space available, FRTs may play 
a role in the ARCADE 2 (and older experiments) excess. 
To that end we have used both an estimate based on 
a differential flux density and a detailed calculation for 
transient sources to calculate the contribution to a radio 
background. 

This analysis has shown that based on what we cur¬ 
rently know about FRBs, they do not contribute at an 
appreciable level to the radio excess measurement. One 
could construct very different source models which could 
give a contribution to (or explanation of) the ERB, and 
such model building is one possible direction to studying 
fast transients and the radio background. 

Overall, this work is meant to serve as a first step 
in a new direction for transient radio sources and the 
unknown cosmic radio background. There is ample op¬ 
portunity here to propose fast radio transient models, 
search and hopefully discover their origin, and perhaps 
link them to the observed, yet unexplained, radio back¬ 
ground. 
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